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ABSTRACT 


A one-dimensional mathematical model for the lithium/thionyl chloride pri- 
mary cell has been developed to investigate methods of improving its performance 
and safety. The model includes many of the components of a typical lithium/thionyl 
chloride cell such as the porous lithium chloride film which forms on the lithium 
anode surface. The governing equations are formulated from fundamental conser- 
vation law r s using porous electrode theory and concentrated solution theory. The 
model is used to predict one-dimensional, time dependent profiles of concentration, 
porosity, current, and potential as well as cell temperature and voltage. When a 
certain discharge rate is required, the model can be used to determine the design 
criteria and operating variables which yield high cell capacities. Model predictions 
can be used to establish operational and design limits within which the thermal 
runaway problem, inherent in these cells, can be avoided. 



INTRODUCTION 


The lithium/thionyl chloride (Li/SOCl 2 ) cell is an attractive primary energy 
source because of its high energy density (1, 2). However, researchers have observed 
that the Li/SOCl 2 cell is a serious safety hazard under certain conditions (2). High 
discharge rates and high temperatures promote thermal runaway which can result 
in the venting of toxic gases and explosion. A mathematical model of this battery 
has been developed to investigate the operational and design characteristics which 
can be adjusted to yield efficient, yet acceptably safe Li/SOCl 2 cells. 

Description of Li/SOCl 2 cell 

The model describes the Li/SOCl 2 cell illustrated in Fig. 1. The four cell 
regions are the lithium foil anode, the lithium chloride (LiCl) film which forms on 
the anode surface, the separator (usually a glass matting), and the porous carbon 
cathode. The anode surface and cathode current collector are the boundaries 
of the model region. The electrolyte consists of lithium tetrachloroaluminate 
(LiAlCLi) in thionyl chloride (S0C1 2 ). The Li/SOCl 2 cell is a complex chemical and 
electrochemical system involving an unknown number of reactions (2). The overall 
reactions which seem most important to include in modeling are: the oxidation of 
lithium at the anode 

Li — > Li + + e~ [1] 

the reduction of S0C1 2 at the cathode 

2S0C1 2 + 4e~ — ► 4Cr + S0 2 + S [2] 

and the formation of the LiCl film on the anode and the precipitation of LiCl on 
the pore surfaces of the cathode during discharge by the following reaction 

Li+ + cr ^ LiCl(s) [3] 
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Reaction [3] is assumed to occur completely and instantaneously, thus, reactions 

[2] and [3] can be combined to yield one overall reaction representing the processes 
occurring within the porous cathode 

4Li + + 4e~ + 2S0C1 2 — 4LiCl + S0 2 + S [4] 

Four species are included in the model: Li + , AlCl^", S0C1 2 , and LiCl. Since the 
chloride ion, Cl - , does not appear in reaction [4] it is not included in the model. 
Later, this species may be included because of its influence on species transport 

(3) . In addition, S0 2 and S are not included in the model because it is assumed 
that they dissolve into the electrolyte (4). These assumptions simplify the model 
development; the more species considered, the more complex and cumbersome the 
transport equations become. The species S and S0 2 may be included at a later 
time because researchers have shown that S0 2 is the major component in the vapor 
phase and that S precipitates with the LiCl in the porous cathode (5). 

Szpak and Venkatsetty (6) state that modeling the Li/SOCl 2 cell is difficult due 
to the complexity of this physical system. They observe that as the cell discharges, 
the temperature and pressure of the cell change, the volume of electrolyte decreases, 
the electrolyte composition varies, and new phases appear. However, much can 
be learned from attempting to describe these phenomena through mathematical 
modeling. 


Models of the Li/SOCl 2 cell 

Previous workers have presented mathematical models of the Li/SOCl 2 cell 
(3, 7-14) as well as mathematical models of similar electrochemical systems (15- 
19). The Li/SOCl 2 cell model developed by Tsaur and Pollard (7) is based on 
conservation of mass and charge and is used to determine the concentration and 
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reaction rate profiles in the cell. Other Li/SOCb cell models (10-14) are based on 
conservation of energy and are used to determine the thermal behavior of the cell. 

Tsaur and Pollard (3, 7, 8) present a one-dimensional model of the Li/SOC^ cell. 
They use porous electrode theory (20, 21) and concentrated solution theory (22-24) 
to develop this model, and then use their model to predict concentration profiles, 
reaction rate distributions in the porous electrode, distributions of LiCl precipitate 
in the pores of the porous electrode, and cell utilization. However, in their model 
development Tsaur and Pollard (3, 7, 8) assume that the reservoir is a well mixed 
solution of constant concentration, which may not be a valid assumption as shown 
by Gu et al. (15). Also, they do not treat the LiCl film as a separate region in their 
model but end their calculations at the LiCl film/separator interface. This work 
shows the importance of including the film region in the modeling as evidenced by 
the large concentration gradients in the film and the effect of film porosity on cell 
performance. 

Several models (10-14) have been presented to describe the thermal behavior 
of Li/SOCl 2 cells. Parnell and Szpak (10) present a thermal model for the 
Li/SOCl 2 cell in one dimension. They consider heat generation due to the 
polarization of the electrodes and internal cell resistance. Their model is used 
to show that as the rate of discharge is increased, the temperature increase is 
proportional to the increase in cell current. Szpak et al. (11) present a thermal 
model, similar to the model by Parnell and Szpak (10), which accounts for 
catastrophic thermal runaway. This model (11) is formulated based on the fact 
that thermal runaway is associated with ignition and burning, via one or more cell 
defects, rather than explosion. The model is used to predict the time and position 
dependent temperature and concentration profiles. The model predictions compare 
favorably to what is physically observed which may indicate that the reactions 
causing thermal runaway are initiated by localized heat sources developed from 
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defective cells and/or cell components. Cho and Halpert (12) have presented a 
simple thermal model for Li/S0Cl2 primary cells. Using experimental data for a 
specific cell discharge to determine a heat generation rate term, they use their model 
to predict the cell temperature as it changes with time. This information can be 
used to determine safe operating limits for a particular cell. Cho and Halpert (13) 
and Cho (14) have taken the modeling one step further by developing a means to 
calculate the heat generation rate term in their model (12) using experimentally 
determined resistances. These thermal models (10-14) are simple models which 
describe the cell using a small set of equations, unknowns, and parameters. These 
models are limited as far as their ability to predict cell behavior as a function of 
various design criteria because they do not include many features of the physical 
system. For example, the models do not include the effects of migration on species 
transport, the effect of the precipitating LiCl in the porous electrode on reaction 
kinetics, and the effects of current and potential on the electrochemical reaction 
rates. These models cannot be used to predict certain key design criteria, such as 
the porosity distribution in the porous cathode. 

Models of similar electrochemical systems, such as the models of the lead acid 
battery (15-17) and the lithium-aluminum/iron sulfide battery (18, 19), were found 
to be useful in the development of the model presented here. These models, like 
Tsaur and Pollard’s model (7), are based on conservation of mass and charge and 
employ concentrated solution theory as well as porous electrode theory. Some of 
the model equations here are similar to those presented in these models. Although 
similar, these models cannot be used to model the Li/SOCU cell because they do 
not include all the various regions in the Li/SOCl 2 cell nor do they include certain 
physical phenomena specific to the Li/SOCl 2 cell, such as the LiCl film formed on 
the lithium anode surface. 
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MODEL DEVELOPMENT 


The governing equations presented next describe conservation of mass and 
current, species transport, and reaction kinetics in each region of the Li/S 0 Cl 2 cell. 
The equations for the porous electrode region are presented first, followed by the 
equations for the reservoir, separator, and LiCl film. The boundary conditions and 
initial conditions necessary to complete the equation set are presented next. The 
overall energy balance used to determine the cell temperature as it changes w T ith 
time is presented last. 


Porous electrode region 

In order to develop a useful and descriptive set of governing equations for the 
porous electrode its essential features must be accounted for without considering 
its random microscopic geometry. The macroscopic theory for porous electrodes, 
introduced by Newman and Tiedemann (20) is used for this purpose. Essentially, the 
porous region is treated as a superposition of two continua, the ionicallv conducting 
electrolyte (solution) phase and the electronically conducting matrix (solid) phase. 
That is, the porous electrode is pictured as one homogeneous region w T hich can 
be described by various distributed quantities. These quantities are continuous in 
space and time and are averaged over a differential volume of the porous region 
which is small relative to the overall dimensions of the electrode yet large relative 
to the pore structure. This volume averaging procedure has been used to model 
other electrochemical systems (15, 25) and is discussed in detail by Trainham (21), 
Whitaker (26), and Dunning (27). Examples of volume averaged quantities are 
the Faradaic current in the solution phase and the electronic current in the matrix 
phase of the porous electrode. 

The current in the porous region must be conserved. Assuming that the 
double layer at the pore wall is extremely thin relative to the pore dimensions, 
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the electrolyte within the pores can be considered electrically neutral and governed 
by the electroneutrality condition which is 


Z{Ci = z + c + + z_c_ = 0 




Then, conservation of charge requires that the charge entering the solution phase, 
i 2 , must be equal to the charge leaving the matrix phase, i 1? which in one dimension 
is expressed as follows 


d% 1 ,z dl2 , x 


= 0 


[e; 


dx dx 

The total current transferred from the matrix phase to the solution phase, di\ tX /dx 
or di 2 ,x/dx, is equal to the net rate of electrochemical reaction per unit volume of 
electrode 

di 2 ,x dij x . , 


= a ik = 


I i 


dx dx 

* 

where j is the local transfer current in A /cm 3 and i* is the rate of electrochemical 
reaction k per unit area of active electrode surface area. 

The material balance equation for each species i in the solution phase is, 

d(eci) 


dt 


= -V • N» + Ri 


(i — , o) 


[ 8 ] 


w T here, 


*• = -£ 


~^F lk 


[9] 


The subscripts +, — , and o refer to Li + , AlCl^, and S0C1 2 , respectively. In Eq. [8], 
e represents the fraction of the porous electrode volume which is void of the matrix 
phase and filled with electrolyte, c,- is the concentration of species i in the solution, 
Nj is the flux of species i averaged over the cross sectional area of the porous 
region, and Ri is the production rate of species i per unit volume of electrode due 
to chemical or electrochemical reaction. Here, reaction [4] is assumed to be the only 
reaction occurring within the porous region. 
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Concentrated solution theory (22-24), in particular the multicomponent dif- 
fusion equation (MDE), is used to develop the equations describing the species 
transport. This theory is used because the soh’ent, SOCI 2 . is involved in an elec- 
trochemical reaction and is not present in excess amounts relative to the solute, 
LiAlCL}. The MDE for a given species can be thought of as a force balance which 
describes how the motion of that species is affected by the motion of all the other 
species in the electrolyte. The MDEs are inverted (22-24) to obtain explicit expres- 
sions for the species fluxes. The resulting flux expressions for a binary electrolyte 
are 

N i = -D eff Vci+ l -^+CiV a (i = +,-) [10] 

z% -t 1 

N c — — D e ff Vc 0 + c 0 v" [11] 

In Eq. [10] and [11] the volume average velocity is defined (28) 

v -=x;N < vi=x:^v < (i=+,-, 0 ) [12] 

i i 

and the effective diffusivity, D e ff , is defined 


D eff = De expe 


[13] 


where expe is chosen to be 1.5; this value is consistent with earlier work (7, 15). Eq. 
[13] relates the electrolyte diffusivity in the porous region, D e ff, to that in the bulk 
solution, D, by accounting for the effects of porosity and tortuosity in the porous 
region (7, 15). 


Solvent and electrolyte material balances: 

Equations [9], [10], [11], and [13] are introduced into Eq. [8] to obtain the 
following material balance expressions for the charged species and for the solvent 


d{ec i) 


= -V 


i ,/■ 

-De expe V Ci + ^+Civ' 

z { F 


V - ' 5 !,fc a • 
- > 7Z l k 

n k F 


dt 


k 


(* = +,-) [14] 



d(tc 0 ) 

dt 


= -V.(- 


De expe Vc 0 


•v - E 


^o 7 k^ 

n k F 




[151 


Equations [14], corresponding to the two unknowns c+ and c_, can be reduced to 
one equation for one unknown by introducing the following relationship 


c = 


v+ 


c_ 

V- 


[16] 


Equation [16] expresses conservation of mass for a completely dissociated salt, here 
LiAlCU, where c is the concentration of the salt (or electrolyte). In Eq. [16], and 
v _ are the number of Li + and AlCl^" ions formed from the complete dissociation of 
one molecule of LiAlCU. respectively. Substituting Eq. [16] into Eq. [14] yields the 
electrolyte material balance 

d(tC) = -V . (-D<'”‘Vc + AfE + cv*') - y [IT] 


at 


Z+V-lF 


^ nv+F 

k 


Note that either relation in Eq. [16] can be substituted into the corresponding 
species balance, Eq. [14], to obtain the same result, Eq. [17j. The cation ( i — -f) 
material balance was chosen here to derive Eq. [17]. 

Solid phases material balance: 

Equation [8] for the solid phases, LiCl here, is 


solid phases * 

|= E E 2 ^-* 


n k F 


[18] 


Equation [18] describes the porosity change of the porous cathode as LiCl precipi- 
tates in the pores. 

Overall material balance: 

The three material balance equations are: Eq. [15] for the solvent, Eq. [17] for 
the electrolyte, and Eq. [18] for the solid species. These equations contain the four 
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unknowns, c, c 0 , e, and v". The unknown c Q can be eliminated from these equations 
via the following independent thermodynamic relationship 

1 =cV e + c 0 V 0 [19] 

where the partial molar volume of the electrolyte, T e , is defined as 

V e = v+ V+ + v- V- [20j 


Equation [19] expresses conservation of volume and can be thought of as a definition 
for the partial molar volumes \\ (29). When Eq. [141 and [15] are premultiplied by 
V+, 1_ , and Vq, respectively, and the three resulting equations are substituted into 
Eq. [18] an overall material balance is obtained. This is true because the resulting 
equation is a linear combination of all the species balances. Equation [19] is used 
to simplify this balance to yield the following overall material balance expression 


Vv 1 


-£ 


n k F 


solid species 

^ 5 i f k 1 ' i “H ( ^ o $ o , k “f~ le^ — — 5 fc 


V i 2 [21] 


In one dimension, Eq. [21] can be integrated to yield an explicit expression for v" 
in terms of i2 lX • Based on the fact that the current is zero when the velocity is zero, 


V x = *2,x = 0 


[ 22 ] 


integration of Eq. [21] yields 


solid species 


V x — ^ jp y y Si^kVi “h o^o.k T 1 e^ — T 1'e^-f.^ — t k 


l 2 ,2 


[23] 
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Ohm’s law for the electrolyte: 

The Faradaic current density, i 2 , is a function of the chemical potential, p e , 
and solution potential, $ 2 (20, 24). This relationship, Ohm’s law for the solution 
phase, for a binary electrolyte is (24) 

«c// r \™4^+ Z+V+ n 4C 0 J 



where 


v = I/_j_ + V- 


[25] 


and 

He = = i/iZT 1 \n(cfa e ) [26] 

In Eq. [24] the effective electrolyte conductivity, /c e // ? is defined analogously to the 
effective electrolyte diffusivity, D e ff. 


Keff = KC 


expe 


[27] 


Equations [19], [26], and [27] can be substituted into Eq. [24] to yield the following 
Ohm’s law expression for the electrolyte 

*2 _ _ V$2 _ vRT { ^ + ,4 SqJCVq \ Vc 

K£ expe 2 F \ra 4 i/+ Z+V+ n 4 (l -cV e )J c 

Ohm ’s law for the matrix phase: 

Ohm’s law for the matrix phase is (15, 20) 



i] = -<7 c// V$i 


[29] 


where 


a eff = <r(l - e°) 


expe 


[30] 
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Equation [30] is used to account for the effect of the matrix phase tortuosity on 
the conductivity of the solid material, a . The initial porosity, e°, is used in Eq. 
[30] because the matrix conductivity is assumed to change insignificantly as LiCl 
precipitates. 

Equation [6] can be used to eliminate ii from Eq. [29]. Integrating Eq. [6] 
with respect to x and applying the boundary conditions 

at x = 0, = I and z 2jZ = 0 [31] 

results in the following relationship 

T ^2,z — I [32 j 

Substituting Eq. [32] into Eq. [29] yields 

I — ii,x — ~ ae ff~Q^~ 1 ^ 3 ] 


Electrode kinetics : 

The electrochemical reaction rates, i*, in Eq. [17] and [18] depend upon 
the local concentrations of the chemical species, the potential driving force for 
reaction, and temperature. Analogous to chemical kinetics, these dependencies 
are not described by fundamental laws. The Butler- Volmer polarization equation 
(15, 30) is used here 


Ik — io,k,ref 





ni 


Co 


href 




exp 


a 


c, k-* 


RT 


-(*1 - ^2 - U hl ref) 


[34] 


The potential driving force for electrochemical reaction k is $1 — $2 — Uk,ref 
where Uk,ref is the equilibrium potential of reaction k measured relative to the 
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reference electrode. The electro-kinetic parameters i 0 ,k,refi (*a,k, and a Ci k are 
the exchange current density, anodic transfer coefficient, and cathodic transfer 
coefficient, respectively. The transfer coefficients for a given reaction sum to the 
total number of electrons involved in that reaction 


a a ,k -r a c .k = «■ k 


[35] 


The rate of reaction k per unit volume of porous electrode is aik where a is the 
available active surface area per unit volume of electrode. The variable a changes 
as the LiCl precipitates on the pore walls, according to reaction [4], and can be 
described by (27, 31) 



where mp is an experimentally determined parameter used to describe the morphol- 
ogy of the precipitate. In writing Eq. [36] it is assumed that all the Li Cl produced 
from reaction [4] precipitates instantaneously and passivates the surface it covers. 
Large values of mp indicate needle-shaped deposits whereas small values represent 
flat deposits (31). 

Substituting Eq. [34] and [36] into Eq. [7] yields the following equation 
expressing conservation of current 


dz 2 , x 
dx 





-h - j n(^3) + 


[37] 


Altogether, six equations and six unknowns are used to describe the porous 
electrode region. The six equations are the electrolyte material balance, Eq. [17]; 
the solid species material balance, Eq. [18]; the overall material balance, Eq. [23]; 
Ohm’s law for the electrolyte, Eq. [28]; Ohm’s law fo- the matrix phase, Eq. [33]; 
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and conservation of current using the Butler- Volmer equation, Eq. [37]. The six 
unknowns are the electrolyte concentration, c; the porosity, e; the volume average 
velocity, u"; the Faradaic current density, z 2 , x ; the potential in the matrix phase, 
$ 1 : and the potential in the solution phase, $ 2 . Simplified forms of Eq. [17], [18], 
[23], [28], [33] and [37] are used to describe the remaining regions of the cell: the 
reservoir, the separator, and the LiCl film. 


Reservoir 

The reservoir, the space between the separator and the positive porous elec- 
trode, contains electrolyte and serves to supply S0C1 2 to the positive electrode as 
it is consumed via reaction [4]. The solid matrix phase does not exist in this re- 
gion, therefore the unknowns e, i ljX , and $ 2 have no meaning. They are treated as 
dummy variables and their values are arbitrarily set as follows 


e = 1 

$i = 0 


^1 ,X — 0 

Using Eq. [40] in Eq. [32] yields 

l2,i = I 


[38] 

[39] 

[40] 


[41] 


Equations [17] and [28] for the porous electrode can be simplified for the 
reservoir using Eq. [38] and [41]. The electrolyte balance becomes 


= nWc - v • (cv - ) 


(42] 


and Ohm’s law for the electrolyte becomes 


k r \n 4 i/. r z+i'-. 


"®o,4 o \ ^ C 

n 4 (l-cl' e )/ C 


[43] 
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Note that effective properties are not used in Eq. [42] and [43] because no matrix 
phase exists. Equation [23] is used to calculate the volume average velocity. In 
summary, the six governing equations for this region are Eq. [23], [38], [39], [41], 
[42], and [43]. 


Separator 

The separator is an inert, porous layer (usually a glass matting material) which 
prevents the negative and positive electrodes from touching and thus shorting the 
cell. The separator porosity (e s ) is constant, hence 

e = e, [44] 

The solid phase of the separator is a non-conducting material; 4>i does not exist 
and is set to zero as in Eq. [39]. All current resides in the electrolyte, hence Eq. 
[41] applies. Using Eq. [41] and [44] to simplify Eq. [IT] and [28], the electrolyte 
balance for the separator is 

Br 

e s - = Der s V 2 c-Y-(cv m ) [45] 

and Ohm’s law for the electrolyte is 

_L = _y<i> 2 + 

K£ e s xs F \n 4 i>+ z+v + n 4 (l - cV e )J c 

Note that effective properties have been used in Eq. [45] and [46] to characterize 
the transport and conductivity in the separator. The exponent exs is analogous 
to expe and its value is set to 1.5. The overall material balance is Eq. [23]. In 
summary, the six governing equations for the separator are Eq. [23], [39], [41], [44], 
[45], and [46]. 


Lithium chloride film 
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The stability of the lithium anode is attributed to the LiCl film which forms 
from the corrosion of the Li by the SOCI2 on open circuit (32-35) 

4Li + 2SOCI2 — ■* 4LiCl + S0 2 + S [47] 

Researchers (32-35) believe that the film passivates the Li and inhibits the mass 
transport of S0C1 2 to the Li surface. Some workers (34, 35) have described the LiCl 
film as consisting of two parts: a relatively thin, non-porous layer of LiCl formed 
on top of the Li anode surface called the solid electrolyte interphase (SEI) and a 
much thicker, porous layer of LiCl formed on top of the SEI called the secondary 
porous layer. These workers have proposed that it is the SEI which passivates the Li 
surface because the secondary layer has a coarse crystalline, porous structure (32) 
which cannot be assumed to have the resistance properties necessary for complete 
passivation of the Li. Peled (35) views the secondary porous layer as a membrane 
which inhibits the flow of ions to and from the SEI. Here, the transition between 
the SEI and secondary layer is assumed to be distinct, and the layers are treated 
separately. The SEI surface is treated here as the Li anode surface having modified 
(passivation) kinetics and the secondary porous LiCl layer is treated here like a 
separator. 

The governing equations for the secondary, porous LiCl film are developed in 
the same manner as those for the separator. The equations are: a constant porosity 

e = t f [48] 


the electrolyte balance, 

Ohm’s law for the electrolyte, 

_i_ - _v$ ( s +,4 , s 0 ' 4 cV 0 \ Vc 

Ke) xf 2 F \n 4 ^ T n 4 (l — c\ e )J c 


[49] 


[50] 
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and Eq. [23], [39], [41]. The value of the exponent exf is set to 1.5. 

The secondary porous film grows as time progresses but its thickness is kept 
constant here. Film growth may occur by reaction [3] where the Cl - ions are a 
product of the cathode reaction. Tsaur and Pollard (7) assume that the rate of 
reaction [47] governs the film growth. However, reaction |47i is the sum of the two 
electrochemical half reactions [1] and [4j and does not occur on discharge because 
the lithium is anodically protected. The secondary porous layer is not of uniform 
thickness over the surface of the lithium (33), its morphology depends strongly on 
the storage temperature (33), and mechanical rupture may occur during discharge 
(32). Because of these complications, the film thickness is kept constant and can 
be thought of as an average value of the thickness which varies in space and time. 
The model will be used to predict the effects of the film thickness and porosity on 
the cell performance. 


Boundary conditions 

The boundaries, from left to right in Fig. 1, are the Li anode surface (taken 
as the SEI surface here), the LiCl film/separator interface, the separator /reservoir 
interface, the reservoir/porous cathode interface, and the cathode current collector 
surface. The governing equations at each boundary are given next. 


SEI surface: 

The SEI is treated as the Li anode surface with modified kinetics. The rate 
of reaction [l] is calculated using Eq. [34] with kinetic parameters representing the 
passivation. Lithium ions must travel from the SEI surface at the same rate at 
which they are produced, thus the flux of this species must equal its reaction rate 


N+ = -u+De e f ef Vc\f 


i 2 r 


z+F 


■ • 

- v+cx = -Zi 


TIjF 


[51] 
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The solvent flux is zero at this boundary 

N o = -De e f xf Vc o \ f + c o x m = 0 [52] 

Using Eq. [19] , Eq. [521 can be written 

N 0 = De e f xf V e Vc\ f + (1 - c \ e ) v* = 0 [53] 

In Eq. [51] - [53], and in equations to follow, the vertical bar and subscript, [ y , 
indicate that the preceding quantity is evaluated at the interface of the region y. 
At the film interface, all the current is transferred from the anode surface to the 
solution phase, hence is zero and Eq. [41] applies. Equation [7] for a flat plate 
electrode is 

'=£'* [54] 

k 

where Eq. [341 is used to calculate t*.. The porosity is treated as a dummy variable 
and it is arbitrarily set equal to the LiCl film porosity, Eq. [48]. At this boundary 
is arbitrarily set to zero, Eq. [39], to establish a reference point from which to 
calculate all other potentials. The absolute values of and $2 have no significance; 
only the value of the difference — $2 is meaningful. The value of $ 2 , given Eq. 
[39 j, is determined from Eq. [54]. In summary, the boundary conditions at the SEI 
boundary are Eq. [39], [41], [48], [51], [53], and [54]. 

LiCl film /separator interface: 

At this boundary the flux of each species i is continuous 



and the velocity is continuous 



[55] 

[56] 
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Substituting Eq. [10], [16], [19], and [56] into Eq. [55] yields the following equation 
expressing continuity of the electrolyte flux 


£ 




[57] 


Equation [46] for the separator and Eq. [50] for the LiCl film are equated, and Eq. 
[57] is used, to yield the following equation which expresses continuity of Faradaic 
current 

= [58] 

As in the film and separator regions, Eq. [39] and [41] apply. The porosity is 
arbitrarily set equal to e s , Eq. [44]. In summary, the boundary conditions at this 
boundary are Eq. [39], [41], [44], [56], [57], and [58]. 


Separator /reservoir interface: 

The boundary conditions here are analogous to those at the LiCl film/ separator 
boundary. These equations are: continuity of flux 


e 


exs 

s 


Vc = Vcl 

£ IT 


continuous velocity 



continuity of Faradaic current 


[59] 

[60] 


ef s y$ 2 | s = V$ 2 | r [61] 

and Eq. [39], [41], and [44]. 

Reservoir/porous electrode interface: 

The boundary conditions at this interface include continuity of flux 

Vc| r = c expe Vc\^ [62] 
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continuity of Faradaic current 


V$ 2 |„ = (“’''V$ 2 | j>e (63] 

and Eq. [41] . This boundary includes the porous electrode surface, hence porosity 
will change as will the volume average velocity as described by Eq. [18] and [23], 
respectively. Equation [33] is applied at this boundary as follows 

V$, I pe = 0 [64] 


Porous electrode current collector surface: 

A Li/SOCl 2 battery may be pictured in one dimension as a series of cells 
sandwiched together. The current collector of the porous cathode separates the 
cathodes of two consecutive cells. With this configuration in mind, the one 
dimensional distribution of the Faradaic current density will be essentially zero 
at the current collector. That is, it is assumed that all of the current has been 
transferred from the solution phase to the matrix phase at the current collector, 
thus 

*2,z = 0 [65] 

Also, the volume average velocity is assumed to be essentially zero at this boundary 

v m x = 0 [ 66 ] 

As at the SEI surface, the flux of the each species must equal its reaction rate. The 
electrolyte flux is zero here 

N+ = -u + De expe Vc\ pe + + v+cv m = 0 [67] 

z +r 

Using Eq. [65] and [66], Eq. [67] can be written 

Vcjpe = 0 [68] 
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Similarly, Eq. [28] becomes 


v$ 2 | pe = 0 [69] 

The transfer current at this boundary is given by 


where j is calculated using Eq. [7], [34], and [36]. The porosity change is calculated 
using Eq. [18]. 


Initial conditions 

To complete the equation set, the initial conditions must be specified for those 
variables that depend explicitly on time, c and e. These initial conditions are 

C — Cinil (all x) [<l] 

e = e° (porous electrode) [72] 


Cell energy balance 


An energy balance is included in the model to determine the cell temperature 
as it changes with time. The energy balance follows from the first law of thermo- 
dynamics 


du dq dw 

dt dt dt 


[73] 


where u is the internal energy of the cell. When constant volume conditions are 
assumed 

du 

dt dt 

The change in the work, dw/di, is zero; no work is done on or by the cell when 
volume is constant. The change in the heat content of cell, dq/dt , is attributed 
to three sources (sinks): heat transferred across the cell boundaries (qht), heat 
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generated by polarization ( q v ), and heat generated by entropy effects ( q s ) (36). 
The heat transfer can be described by 


qht = -h 0 (T - T a ) [75] 

where h a is the heat transfer coefficient characterizing the heat flow from the cell to 
the surroundings. Here, the temperature of the surroundings is assumed to remain 
constant at Ta- The remaining two sources of heat are given by (36) 

q p + q s = I^E oc - E - T-^j = I{Etn ~ E) [76] 

E tn in Eq. [761, termed the thermoneutral potential, is the theoretical open circuit 
potential of the cell at absolute zero. Substituting Eq. [74], [75], and [76] into Eq. 
[73] yields the overall energy balance 

C ^ = —h 0 (T - T a ) + I(E tn - E) [77] 

Note that each term in Eq. [77] is based on the total active surface area in the cell. 
The initial condition for T here is 

T = T A (all x) [78] 
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MODEL PARAMETERS 


The parameter values input to the model are listed in Table I. The partial 
molar volumes and transference numbers were assumed to be constant throughout 
the discharge at the values listed in Table I. The kinetic parameters i 0) i iT . e /, oc a ,i, 
a c , i, and a°i 0 , 4 ,re/ were chosen so that the model predictions compared favorably 
with experimental data (32, 37-41). A temperature dependent expression for the 
diffusion coefficient and a temperature, concentration dependent expression for the 
conductivity are given in the appendix. 

The independent adjustable parameters in this model development are the 
porous electrode thickness ( S pe ), the initial porosity of the porous electrode (e°), 
the thickness of the electrolyte reservoir ( S r ), the thickness of the separator (S s ), 
the porosity of the separator (e s ), the LiCl film thickness ( Sj ), the LiCl film 
porosity (e/), the initial elect rolyte concentration (c; n j t ), the heat transfer coefficient 
( h a ), the current density (7), and the temperature of the surroundings (Tjl). One 
objective of this work is to use the model to determine which of these independent 
adjustable parameters affects cell performance the most. Table II lists typical values 
for these independent adjustable parameters (5). Unless otherwise noted, these base 
case values are used throughout the work presented here. 
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METHOD OF SOLUTION 


The system of coupled, nonlinear, ordinary differential equations describing the 
Li/SOCU cell, shown in Fig. 1, is solved numerically because an analytical solution 
is not possible. The spatial derivatives are approximated using finite differences 
accurate to O(Ax) 2 where Ax is the largest step size in the x direction. Implicit 
stepping, accurate to 0(At ), is used for the time derivatives. A Newton-Raphson 
algorithm using deBoor's banded matrix solver (42-44) is used to solve the system 
of finite difference equations. The procedure is iterative and requires initial guesses 
of the unknowns. Step sizes in both space and time w'ere reduced until predicted 
cell voltages were reproducible to three significant digits. 
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RESULTS AND DISCUSSION 


Figure 2 shows predicted electrolyte concentration profiles end how they 
change with time. Note that the concentration profiles are relatively flat in the 
separator and reservoir regions indicating that diffusion, migration, and convection 
of species are essentially uninhibited in these regions. Convection contributes least 
to the species transport; v" is O(10~ 5 ) and the convection term is an order of 
magnitude less than the diffusion and migration terms in Eq. [17], [42], [45], 
and [49]. The large concentration gradients in the film region are due the small 
porosity of this region w’hich reduces significantly the effective diffusivity. Near 
the beginning of discharge transport in the porous electrode region is essentially 
uninhibited. As the discharge proceeds, LiCl begins to fill the pores of the electrode 
and transport becomes more difficult as evidenced by the steeper concentration 
gradients at later times. The electrolyte concentration increases throughout the 
cell as time progresses and the solvent concentration decreases according to Eq. 
[19]. These profiles show that throughout most of the discharge cell performance is 
dictated by the rate of reaction [4] at the cathode. Tsaur and Pollard (7) predict 
diffusion control as evidenced by the large concentration gradients in the separator 
of their simulated cell; they assume a constant reservoir concentration a priori. 
This difference, reaction control here versus diffusion control (7), could be due to 
the different governing equations and boundary conditions as well as the differing 
input parameters used in their model and the present model. From Fig. 2 it 
can be seen that the change in the concentration across the reservoir is about half 
what it is in the separator and, under conditions of diffusion control, this change 
could effect results significantly. Higher discharge rates, lower temperatures, faster 
porous electrode reaction kinetics (i.e., larger exchange current densities for the 
thionvl chloride reduction on the carbon electrode), lower initial concentrations, 
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and smaller porosities all contribute to changing the reaction controlled conditions 
observed here to diffusion controlled conditions. 

Figures 3-5 illustrate the space-time behavior of the porous electrode porosity, 
Faradaic current density, and overpotential, respectively. These profiles show that 
the reaction rate distribution changes very little as time progresses. The reaction 
rate will be greatest at the front of the porous electrode because the matrix 
conductivity (<r) is much larger than the electrolyte conductivity («■) (7). That 
is, at the reservoir /'porous electrode interface all current resides in the solution 
phase but upon entering the porous electrode the current will be transferred 
rapidly to the matrix phase because the matrix phase offers less resistance to 
current flow. However as LiCl precipitates, fills the pores as shown in Fig. 3, and 
decreases the available active surface area, according to Eq. [36], the reaction rate 
distribution becomes more uniform. This change can be seen in Fig. 4 where the 
i 2 , x distribution flattens out; the reaction rate is the slope if the % 2 , x versus x curve. 
Figure 5 shows that the overpotential, the driving force for electrochemical reaction, 
is highest at the front of the electrode as expected. The overpotential increases with 
time indicating that the electrochemical reaction must proceed at a faster rate to 
compensate for decreasing available active surface area. 

As the end of discharge is approached the front portion of the electrode becomes 
clogged with LiCl precipitate (see Fig. 3). Lithium ions are unable to cross this 
barrier hence the electrolyte concentration increases rapidly in the film, separator, 
and reservoir as it is continually produced at the anode. At this point the electrolyte 
contained in the pores of the porous electrode supplies the lithium ions consumed 
in reduction reaction [4]. The electrolyte concentration drops in the vicinity of 
the highest reaction rate (the front, of the electrode) and the electrolyte supply 
contained in the back portion of the electrode must diffuse to this reaction plane. 
This behavior can be seen from Fig. 4 where the current density distribution at t = 
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2160 seconds falls back from the flatter profiles exhibit at earlier times. Eventually, 
the front of the electrode becomes depleted of electrolyte and is unable to support 
further reaction, hence current flow, and the cell fails. This failure is indicated by 
the model when convergence cannot be achieved. 

The model can be used to demonstrate the effects of the independent design 
criteria on cell performance. The least influential parameters, for the reaction 
controlled conditions here, were found to be S r and S s . For example, decreasing 
either S r or S s by an order of magnitude resulted in changes in the discharge time 
and cell voltage of less than 1%. The most influential parameters on cell performance 
were found to be S pe , e°, and c; n ^. The parameters characterizing the LiCl film, e/ 
and Sf , show small to moderate effects on cell performance. 

As shown in Fig. 6, beyond a porous electrode thickness of approximately 0.1 
cm there is little gain in cell lifetime. Apparently, for thicknesses beyond 0.1 cm the 
advantages of increased active surface area and electrolyte storage are diminished 
by the increase in cell resistance and hence decrease in cell lifetime. It is difficult to 
determine the effect of the porous electrode porosity on cell performance because 
when porosity is changed the active surface area, a, is also changed. Keeping a 
constant here, the cell lifetime is decreased by 49% and 86% for e° decreased to 0.6 
and 0.35, respectively. This decrease in lifetime occurs simply because there exists 
less void volume for the precipitating LiCl to fill. Changes in S pe and e° have very 
little effect on the operating voltage and temperature; the voltage and temperature 
varied by less than 1% for the cases studied here. 

Figure 7 shows the predicted effect of initial electrolyte concentration on 
cell lifetime. Caution is needed when comparing cell performance for different 
concentrations because the average operating voltage deer, ases as initial electrolyte 
concentration decreases. For example, the average operating voltage for Ci n n = 
0.25 M was about 8% lower than that for the base case. This concentration could 
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not be included in Fig. 7 because the predicted voltage transient fell below the 
cutoff voltage (2.9 V) chosen to determine cell lifetime. Beyond a concentration of 
approximately 2 M there is little increase in cell lifetime or average cell voltage. At 
concentrations below 1 M the overpotential required to drive the electrochemical 
reaction at a fast enough rate to support the current draw is necessarily high, 
especially at the front of the porous electrode. This acceleration of reaction rate 
causes the pores to fill up quickly with LiCl, and the cell fails rapidly as shown in 
Fig. 7. Also, for initial concentrations of about 0.5 M and below the shape of the 
cell voltage transient changes; cell voltage initially increases, reaches a maximum, 
and then decreases. This initial increase in voltage is attributed to the rather large 
initial increase in conductivity as the concentration of lithium ions increases as 
they are produced at the anode. Later the voltage decreases due to the decreasing 
porosity in the porous electrode and the resulting increase in resistivity. 

The porosity of the LiCl film which forms on the lithium anode has a small effect 
on cell performance. Porosity of this film varies depending on storage temperature; 
higher storage temperatures tend to produce larger LiCl crystals and more compact 
films (33). Several observations can be made here from the model predictions. Film 
porosities beyond 0.1 seem to have little to no effect on cell lifetime, voltage, or 
temperature. Below a porosity of 0.1, cell voltage begins to decrease due to the 
greater hindrance of species transport in the film. For example, for a porosity of 
0.05, the average cell voltage and cell lifetime are reduced by 3% as compared to 
the base case. When porosity is extremely low (e.g., 0.005) the model predicts 
that the cell will not operate for the 30 mA/cm 2 drain considered here. For a film 
porosity of 0.025, the model predicts higher cell temperatures (about 10°C higher 
on the average) and an average cell voltage about 11% lower than that for the base 
because of the increased cell resistance. These results show 7 that it is advantageous 
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to maintain low temperature storage (room temperature) so that the film does not 
become too compact. 

The film thickness has little effect on cell performance for the conditions studied 
here. Changes in cell lifetime of less than 0.1% are predicted when 5/ is decreased 
to 0.5 microns or increased to 40 microns. However, depending on the film porosit} 7 , 
the average cell voltage can be decreased significantly when the film is thick. For 
the 40 micron film thickness, the average cell voltage was decreased by about 4% 
as compared to the base case. 

The effect of discharge rate on cell voltage and temperature is shown in 
Fig. 8 and 9. Figure 8 shows the expected behavior of increased average cell 
voltage, lifetime, and capacity with decreased discharge rate. These curves compare 
favorably with experimental voltage-time curves (37, 38, 41). It is difficult to make 
detailed comparisons of model predictions with the experimental data reported in 
the literature due to the lack of cell specifications given (e.g., porous electrode 
porosity and thickness, and separator porosity and thickness). The temperature- 
time curves for various discharge rates are shown in Fig. 9. The initial rise in 

temperature is due to the polarization of the electrodes and entropy effects as 

represented by the terms q p and q s in Eq. [76]. As the temperature of the cell 
increases the rate of heat transfer to the surroundings increases until an equilibrium 
is almost reached where the heat transfer almost equals the heat generation. These 
equilibriums are represented by the plateaus in the curves in Fig. 9. At higher 
discharge rates these plateau become less distinct because the heat transfer cannot 
keep pace with the greater heat generation produced by the more severe electrode 
polarizations at higher discharge rates. At the end of discharge localized reactant 
depletion causes severe polarization and the cell temperature is driven upward at 

a faster rate as shown by the tail ends of the curves in Fig. 9. This surge in 

cell temperature at the end of discharge could lead to unsafe behavior by rapidly 
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increasing the internal pressure or possibly triggering exothermic runaway reactions 
(11). To avoid this hazard the cell could be operated below certain discharge rates 
and ambient temperature limits or it could be used up to a limited time well before 
reaching the end of discharge. These temperature profiles exhibit the same general 
pattern of experimental temperature-time curves (39). Improved heat transfer paths 
could be investigated using the model by changing the energy balance parameters, 
C and h 0 , or modifying the energy balance, Eq. [77]. 

Model predictions are only as good as the input to the model. The results 
presented here are in good agreement with experimental data reported in the 
literature perhaps because of the detailed theoretical development and the trial 
and error perturbation of some of the kinetic parameters mentioned earlier. The 
kinetic parameters and physical properties used in the present analysis may change 
significantly for different applications which, in turn, could significantly effect 
model predictions. The versatility to change model inputs to simulate different 
applications is an important aspect of this model. 

The kinetic parameters, f 0j j, re /, a°i 0t 4 jre f. a a> j, a C; i, a aj4 , a c ,4, and the 
morphology parameter mp which is used to calculate a, are the least known 
because these parameters cannot be determined easily by experiment. Also these 
parameters are not constant values; the kinetic parameters for the SOCI 2 reduction 
are specific to the type of carbons used for the porous electrode, for example. 
The morphology parameter will change depending on the discharge rate and cell 
temperature. Experimental data along with complete cell specifications (including 
porosities) and physical property data (such as diffusivity and conductivity) could be 
used in conjunction with a parameter estimation technique to estimate these kinetic 
parameters. A detailed sensitivity analysis of the model predictions to changes in 
these kinetic parameters is beyond the scope of this work though a few items are 
of note. The oxidation exchange current density, i 0 ,i,ref, can significantly effect 
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cell voltage and lifetime whereas the reduction exchange current density, a°i 0 ,4,refi 
seems to influence cell lifetime more than cell voltage. For example, increasing 
a °io,4,ref by an order of magnitude results in a 15% reduction in cell life and a less 
than 1% increase in cell voltage. 

Some of the physical properties used here were estimated from the available 
data in the literature (6, 7, 37, 45-47). An approximate expression for diffusivity, 
as described in the appendix, was used and constant transference numbers were 
assumed. The sensitivity of model predictions to small changes in these properties 
was briefly investigated. If D is halved, cell lifetime decreases by only 6%. If t" 
is decreased to 0.25 there is an 18% decrease in cell lifetime due to the resulting 
decrease in species transport, Eq. [10], and the increase in IR drop across the cell 
as indicated by Ohm’s law for the electrolyte. Changes in D and f" have negligible 
effect on cell voltage and temperature. Similarly, changes in conductivity will 
directly effect the IR drop across the cell as dictated by the Ohm's law expression 
for the electrolyte. Decreasing conductivity will increase the IR drop and decrease 
cell voltage and lifetime. 
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CONCLUSIONS AND RECOMMENDATIONS 


A one-dimensional model of a complete Li/SOCl 2 cell has been developed and 
used to demonstrate both design and operational changes that can be made to 
increase cell capacity and life. An unending number of cases could be studied using 
the model to determine optimal and safe designs for various cell specifications and 
discharge rates. This ability to investigate different designs simply by changing 
model inputs is where the value of the model lies. In addition, the model can be 
used along vcith a parameter estimation technique to estimate kinetic parameters 
for various carbon materials and electrolyte chemistries. 

Improvements in the model would include accounting for the growth of the Li Cl 
film as the cell discharges, the addition of the Cl", SO 2 , and S species to determine 
their effects on cell performance, changing the energy balance to account for more 
complex heat transfer paths, and the modelling of a stack of cells representing the 
cross section of an entire Li/SOCU battery. Also porous electrode swelling has been 
shown to be an important aspect of cell performance during high rate discharge (48) 
and should be added to the model. 
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APPENDIX: PHYSICAL PROPERTIES 


The conductivity of the electrolyte, k, is a function of concentration and 
temperature. Szpak and Venkatasetty (2) suggest that this functionality should be 
included in the battery modelling because of the large changes in the conductivity 
as the battery discharges. An empirical expression, similar in form to other 
conductivity expressions (7, 15, 16), describing this functionality has been developed 
from the available data (6, 37, 45-47) 

f ft\ c exp (ft 2 c + ft 3 c 2 ) exp (ft 4 /T) (c < 1.8 M) [7Ql 

\ /3 5 exp (ft 4 /T) (c > 1.8 M) ‘ /yj 

where ft\ = 135.45 £l~ 1 cm~ 1 , ft 2 = 399.09 cm 3 /mol, ft 3 = -2.5055 x 10~ 5 cm 6 /mol 2 , 

ft 4 — -711.69 K, and ft h — 0.22244 The data for LLAICI4 mixtures 

(equimolar mixtures of Li Cl and AICI3) are for one temperature only, 298.15 K 

(37, 45). The concentration range of this data is 0.000313 M to 2.0 M. The 

dependence of conductivity on temperature and concentration for concentrations 

above 2.0 M was estimated from conductivity data for a similar electrolyte, 

SOCI2 containing S0 2 , S, and equimolar amounts of LiCl and AICI3 (46-47). 

Equation [79] predicts continuous conductivities and derivatives (with respect to 

concentration and temperature) unlike a conductivity expression presented earlier 

( 7 ). 

No data exists to perform a similar analysis for the diffusion coefficient. 
The diffusivity was assumed to be a function of temperature and independent of 
concentration. The dependence of diffusivity on temperature was assumed to be 
identical to the dependence of conductivity on temperature as was done in previous 
work (7). The expression used here is 

D = ft 6 exp (fa/T) [80] 

where ft 6 — 1.0 x 10 -4 cm 2 /s. 
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LIST OF SYMBOLS 


a 
a° 
a 0 
c 

C-i } ref 

C 

D 

D ef f 

E 

Eoc 

E tn 

exf 

expe 

exs 

f 

F 

h 0 

ik 

io,k,ref 

ii 

h 

I 

j 

K it i 

mp 

nk 

N* 

Pi, k 
Qi,!* 

1 

qhi 

Qp 

q s 

R 

Ri 

Si 

$i,k 


specific active surface area of the porous material. cm ~ 1 
initial value of a, cm ~ 1 

property expressing secondary reference state of electrolyte, cm^ /mol 

concentration of the electrolyte, mol /cm 3 

concentration of species i, mol/ cm 3 

reference concentration of species z, mol /cm 3 

heat capacity of the cell, J /cm 2 K 

diffusion coefficient of the binary electrolyte, cm 2 j s 

effective diffusion coefficient of electrolyte, cm 2 ! s 

cell voltage, V 

open-circuit voltage of the cell, V 
thermoneutral potential of t e cell, V 

exponent for calculating effective properties in the LiCl film 
exponent for calculating effective properties in the porous 
electrode region 

exponent for calculating effective properties in the separator 

mean activity coefficient of electrolyte 

Faraday : s constant, 96487 C /mol of electrons 

heat transfer coefficient, W/cm 2 K 

reaction rate of electrochemical reaction fc, A/cm 2 

exchange current density of reaction k at c re f, Afcm~ 

superficial current density in the matrix phase, A /cm 2 

superficial current density in the solution phase, A /cm 2 

total cell current density, A /cm 2 

current transferred between phases, A/ cm 3 

friction coefficient for interaction between species i and L Js/cm 5 
morphology parameter 

number of electrons transferred in the electrochemical reaction k 

superficial flux of species i, mol/ cm 2 s 

anodic reaction order of species i in reaction k 

cathodic reaction order of species i in reaction k 

net heat in cell, J /cm 3 

energy change due to heat transfer between cell and 
surroundings, W /cm 2 

heat generated due to cell polarization, W/cm 2 
heat generated due to entropy effects, W/cm 2 
universal gas constant, 8.3143 J /molK 
production rate of species i due to reaction, mol/ cm 3 s 
thickness of region /, cm 

stoichiometric coefficient of species i in reaction k 
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t 

*: 

T 

T A 


u 

t fc.re/ 


Vi 

v" 

Vi 

w 

X 

Z{ 


time, s 

transference number of species i relative to v" 
cell temperature, K 
ambient temperature. K 
internal energy of the cell, J / cm 3 

Potential of reaction k relative to the reference electrode, V 

superficial velocity of species i, cm/s 

superficial volume average velocity, cm/s 

partial molar volume of species i, cm 2 /mol 

net work done by or on the cell, J /cm 3 

normal direction from SEI surface to porous electrode, cm 

charge number of species i 


Greek Symbols 

a ay k transfer coefficient in the anodic direction of reaction k 

a Ci k transfer coefficient in the cathodic direction of reaction k 

e porosity or void volume fraction of the porous electrode 

e° initial porosity of the porous electrode 

k conductivity of the solution, fl~ 1 cm~ 1 

K e ff effective electrolyte conductivity, n -3 cm -3 

fii chemical potential of species i, J /mol 

v+ number of cations produce by dissociation of electrolyte 

v- number of anions produce by dissociation of electrolyte 

v total number of ions into which one molecule of the salt dissociates 

a conductivity of the solid, ft /cm 

<r e ff effective conductivity of the matrix phase, ft /cm 

<J>j potential in the matrix phase, V 

potential in the electrolyte, V 


Subscripts 

e electrolyte 

/ Li Cl film 

i species i 

init initial value 

k reaction k 

0 solvent 

pe porous electrode 

r reservoir 

s separator 

-f~ cation 

— anion 

1 solid matrix phase 

2 electrolyte phase „ 
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Table I. Parameter values used in model 




Physical properties 



Parameter 

Value 

Ref. 

Parameter 

Value 

Ref. 


i 


z/_ 

1 


z + 

i 


Z __ 

-1 


<■ 

0.5 


t m _ 

0.5 


V e (cm 3 /mol) 

77.97 

7, 38 

V 0 (cm 3 ! mol) 

72.63 

46 

liiCl (cm 3 j mol) 

20.5 

49 

a (fl~ 1 cm~ 1 ) 

45.5 

50 

D (cm 2 / s') 

t 


k (fi -1 ™ -1 ) 

t 


E tn (V) 

-3.723 

36 

dE/dT (V/K) 

2.28x10 

- 4 36 

C (J/cm 2 K) 

0.2 

5, 7, 14 






Kinetic parameters 



Li oxidation (reaction 1) 


SOCI 2 reduction (reaction 4) 

Parameter 

Value 

Ref. 

Parameter 

Value 

Ref. 

*+,i 

-1 


S-,4 

-2 





o,4 

-1 





<SLiCl,4 

2 


9+,i 

1 


9+ ,4 

1 





Qo, 4 

0.5 


io,l, ref (A/cm 2 ) 

5.0xl0~ 4 


a°«o,4,re/ (A/ cm 3 ) 

1.0 





mp 

0.05 

31 

t'w (V) 

E tn + T % f 36 

Ua .re/ (V) 

0.0 


a Q , 1 

0.25 


^a.4 

1 


Q c ,i 

0.75 


«c,4 

1 


n i 

1 


714 

2 


f see appendix 



Table II. Base case values of the independent adjustable parameters 


Parameter 

Value 

Parameter 

Value 

Spe (cm) 

0.0335 

e° 

0.85 

S r (cm) 

0.01 



S s (cm) 

0.0127 


0.7 

Sf (cm) 

0.001 

e f 

0.1 

Cinit ( M ) 

1.5 

h a (J / cm? Ks) 

6.0x10 

I (A /cm 2 ) 

0.03 

Ta (A') 

298.15 
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